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Popular Summary 

The global hydrological cycle is central to the Earth’s climate system, with rainfall and the 
physics of its formation acting as the key links in the cycle. Two-thirds of global rainfall occurs 
in the Tropics. Associated with this rainfall is a vast amount of heat, which is known as latent 
heat. It arises mainly due to the phase change of water vapor condensing into liquid droplets; 
three-fourths of the total heat energy available to the Earth's atmosphere comes from tropical 
rainfall. In addition, fresh water provided by tropical rainfall and its variability exerts a large 
impact upon the structure and motions of the upper ocean layer. 

An improved convective -stratiform heating (CSH) algorithm has been developed to obtain the 
3D structure of cloud heating over the Tropics based on two sources of information: 1) rainfall 
information, namely its amount and the fraction due to light rain intensity, observed directly 
from the Precipitation Radar (PR) on board the TRMM satellite and 2) synthetic cloud physics 
information obtained from cloud-resolving model (CRM) simulations of cloud systems. The 
cloud simulations provide details on cloud processes, specifically latent heating, eddy heat flux 
convergence and radiative heating/cooling, that are not directly observable by satellite. The new 
CSH algorithm-derived heating has a noticeably different heating structure over both ocean and 
land regions compared to the previous CSH algorithm. One of the major differences between 
new and old algorithms is that the level of maximum cloud heating occurs 1 to 1.5 km lower in 
the atmosphere in the new algorithm. This can effect the structure of the implied air currents 
associated with the general circulation of the atmosphere in the Tropics. 

The new CSH algorithm will be used provide retrieved heating data to other heating algorithms 
to supplement their performance. 
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Abstract 

The relationship between the surface rainfall intensity and its associated stratiform amount is 
investigated by examining the Tropical Rainfall Measurement Mission (TRMM) 
Precipitation Radar (PR) product and the Cloud Resolving Model (CRM) simulated 
precipitation data. The results showed that there are more stratiform rain amount for light 
rain for both ocean and land region, but more convective rain amount for heavy rainfall. The 
differences between oceanic and land region is that more convective rain for precipitation 
rate over 100 mm/day for the land region than its oceanic counterpart. 

An improved convective -stratiform heating (CSH) algorithm is developed based on 
the PR rainfall product (rainfall intensity and stratiform percentage) and CRM simulated 
latent heating, eddy convergence/divergence and radiative heating/cooling. The new CSH 
algorithm derived heating has very different heating structure for the ocean and land regions 
compared to previous CSH algorithm derived. Another major difference between new and 
previous CSH algorithm is that the level of maximum heating is 1 to 1.5 km lower for the 
new than old CSH algorithm. 
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1. Introduction 

Latent heat is dominated by phase changes between water vapor and small liquid or frozen 
cloud-sized particles. It consists of condensation of cloud droplets, evaporation of cloud 
droplets and raindrops, freezing of raindrops, melting of snow and graupel/hail, and the 
deposition and sublimation of ice particles. In addition, eddy heat flux convergence from 
cloud motions can also result in distributing the heating or cooling vertically and 
horizontally. Latent heating cannot be measured directly with current techniques. These 
processes are not directly detectable with current remote sensing or in situ instruments, 
which explain why the retrieval schemes to be discussed depend heavily on some type of 
cloud-resolving model (CRM). But, it can be derived indirectly by measuring the vertical 
profiles of temperature and the 3D wind fields from extensive rawinsonde networks through 
a residual method (called a diagnostic budget; Yanai et ol. 1973; House 1982, 1997; Johnson 
1984, Zhang and Lin 1997, Schumacher et al. 2006). 

The launch of the Tropical Rainfall Measuring Mission (TRMM) satellite, a joint 
U.S. -Japan project, in November of 1997 made it possible for quantitative measurements of 
tropical rainfall to be obtained on a continuous basis over the entire global tropics. TRMM 
provides a much needed and accurate measurement of rainfall as well as an estimate of the 
four-dimensional structure of latent (diabatic) heating over the global tropics. In the last few 
years, it has been established that standard products of LH from TRMM measurements 
would be a valuable resource for scientific research and applications (see a review by Tao et 
al. 2006). Such products would enable new insights and investigations concerning the 
complexities of convection system life cycles, the diabatic heating controls and feedbacks 



related to meso-synoptic circulations and then forecasting, the relationship of tropical 
patterns of LH to the global circulation and climate, and strategies for improving cloud 
parameterizations in environmental prediction models. 

Five different TRMM LH algorithms designed for applications with satellite- 
estimated surface rain rate and precipitation profile inputs have been developed, 
intercompared, validated, and applied in the past decade (see Tao et al. 2006; 2007). They 
are the: (1) Goddard Convective-Stratiform Heating (CSH) algorithm, (2) Hydrometeor 
Heating (HH) algorithm, (3) Goddard Profiling Heating (GPROF Heating) algorithm, (4) 
Spectral Latent Heating (SLH) algorithm, and (5) Precipitation Radar Heating (PRH) 
algorithm. The CSH algorithm only requires information on the surface precipitation rates, 
amount of stratiform rain, and the type and location of the observed cloud systems (Tao et al. 
1993). A lookup table, however, is used containing stored convective and stratiform latent 
heating profiles, normalized by total surface rain rates, for various types of cloud systems in 
different geographic locations. These profiles are mostly obtained from CRM (GCE model) 
simulations. In GPROF, CRM-simulated hydrometeor and associated latent heating vertical 
profiles that have radiative characteristics consistent with a given set of multi-spectral 
microwave radiometric observations are composited to create (retrieve) a best estimate of the 
observed profiles (Olson et al. 1999, Grecu et al. 2009). The HH and algorithm 1 estimates 
latent heating profiles as a function of the vertical derivative of the retrieved hydrometeor 
profiles (Yang and Smith 1999). The derivation and evaluation of the HH algorithm is also 
based on CRM simulations. It requires information about the vertical profiles of cloud- and 

1 Tao et al. (1990) were the first to put forth a latent heating algorithm that was termed a 

hydrometeor/lieating algorithm. Tao et al. (1993) then improved the performance of the HH 
algorithm by including surface rain rates. 
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precipitation-sized water and ice particles, all of which can be obtained from the TMI 
profiler retrievals (Smith et al. 1992, 1994). The terminal (fall) velocities of the large cloud 
(precipitating) particles (rain, snow and graupel/hail) are also required for the HH algorithm. 
Cloud-scale velocity is needed in Yang and Smith (1999a,b) and is obtained by applying a 
regression method to a CRM-simulated database. The SLH algorithm (Shige et al. 2004, 
2006, 2009) is also based on GCE model results. It uses PR information (melting layer, 
precipitation top height, rain rate and type) to select the heating profiles in a look-up table. 
The PRH algorithm (Satoh and Noda 2001, Katsumata et al. 2009) also uses PR information 
but without using any heating profiles simulated by CRMs. However, it needs to estimate the 
cloud drafts and (standard) thermodynamic structures associated with cloud systems. An 
iteration calculation is applied to match the relationship between rainfall and latent heating. 
The strengths and weaknesses of each algorithm were discussed in Tao et al. (2006). Table 1 
lists the data required, heating product of each algorithm and key references of these 
algorithms. 

In this paper, both TRMM PR and CRM simulated results will be used to examine 
and quantify the convective and stratiform percentage and their relationship with surface 
rainfall and latent heating for convective systems that developed in different geographic 
locations. The results from the GCE model simulations will be used to calculate the each 
components of atmospheric budget (latent heat, eddy transport and radiation). These 
simulated results will be used for the improvement of a new CSH algorithm. In Section 2, 
the model and model set-up used, and cases are described. The results are shown and 
discussed in Section 3. In section 4, the performance of new CSH algorithm will be 
presented. Finally, the major results and future work are summarized in Section 5. 
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2. Numerical Modeling and Data 

2.1 Goddard Cumulus Ensemble (GCE) Model 

The Goddard Cumulus Ensemble (GCE) model, a Cloud Resolving Model (CRM), is used to 
simulate clouds/cloud systems and their associated heating budget. It is non- hydrostatic. It 
takes account of both the absorption and scattering for solar radiation and both the emission 
and absorption for infrared radiation. Its cloud-radiation interaction has been assessed (Tao et 
al. 1996; 2003). Subgrid- scale (turbulent) processes in the model are parameterized using a 
scheme based on Klemp and Wilhelmson (1978) and Soong and Ogura (1980). The effects 
of both dry and moist processes on the generation of subgrid-scale kinetic energy have been 
incorporated. The sedimentation of cloud ice (Starr and Cox 1985) is included to better 
model clouds in the upper troposphere. All scalar variables (temperature, water vapor, and all 
hydrometeors) are calculated with a positive definite advection scheme (Smolarkiewicz and 
Grabowski 1990). Results from the positive definite advection scheme are in better 
agreement with observations for tropical cloud systems (Johnson et al. 2002). 

The model has five prognostic hydrometeor variables: the mixing ratios of cloud 
water, rainwater, cloud ice, snow and graupel. It has recently modified by Lang et al. (2007) 
to reduce un-realistic large precipitating ice particles and by Zeng et al. (2008, 2009) to 
introduce the ice nuclei concentration into the parameterization of the Bergeron process as an 
input factor. The development and main features of the GCE model were also published in 
Tao and Simpson (1993) and Tao et al. (2003). A review on the application of the GCE 



model to better understand precipitation processes can be found in Simpson and Tao (1993) 
and Tao (2002). 
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In the GCE model, each grid point is designated as either a cloudy or clear area for 
each integration time, depending upon whether the sum of the cloud water and ice mixing 

• Q C 

ratios are larger than 10" (or 10"' depending on the type of clouds/cloud systems) g/kg at 
each grid point (usually an indicator of 100% relative humidity). In the cloudy area, each 
grid point can also be designated as either being in the active or 11011 -active updraft or 
downdraft region (Tao et a /. 1987). I 11 addition, the cloud characteristics can be divided into 
convective and stratiform components (Tao et al. 1993; Lang et al. 2003). In short, 
convective regions include those with large vertical velocities (exceeding 3-5 m s" 1 ) and/or 
large surface precipitation rates. The stratiform region is separated into regions with or 
without surface rainfall. 

The model in the present paper has the same structure as that in previous studies (e.g., 
Soong and Tao 1980; Tao and Soong 1986; Johnson et al. 2002; Tao et al. 2003; 2004; Zeng 
et al. 2007 and many others) wherein clouds are simulated under prescribed large-scale 
forcing. The default numerical experiment is two-dimensional, using a 1 km horizontal 
resolution and vertical resolution that ranges from 42.5 m at the bottom to 1 km at the model 
top, which is at 22.5 km. The model uses a time step of 6 seconds and 512x41 grid points for 
integration. Please see more details in Zeng et al. (2008, 2009). 


2.2 Data 
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(a) Oceanic Convective Systems (GATE, TOGA COARE and SCSMEX) 

The South China Sea Monsoon Experiment (SCSMEX) was conducted in May-June 1998. 
Two major convective events, prior to and during monsoon onset (May 18-26, 1998) and 
post monsoon (June 2-11, 1998) are observed. The SCSMEX forcing data were obtained 
from a variational analysis approach (Zhang and Lin 1997; Zhang et al. 2001). Its numerical 
simulation starts at 0600 UTC 6 May 1998 and lasts for 44 days. The Tropical Oceans 
Global Atmosphere (TOGA) - Coupled Ocean Atmosphere Response Experiment (COARE) 
was conducted in 1992 and 1993 over the Central Pacific. The most intense convection 
during TOGA COARE occurred in middle and late December 1992, prior to the peak of 
westerly wind bursts around 1 January 1993. Several major convective events occurred 
around 11-16 and 20-25 December 1992, mainly due to low-level large-scale convergence of 
easterlies and westerliey (Lin and Johnson 1996). For TOGA COARE, the large-scale 
forcing used in the GCE model was derived from Intensive Flux Array (UFA) sounding 
networks (Lin and Johnson 1996). GARP (Global Atmospheric Research Program) Atlantic 
Tropical Experiment (GATE) was conducted in 1974 over East Atlantic. Cloud systems 
( 11011 -squall clusters, a squall line, and scattered convection) for the period 1-8 September 
1974 phase III of GATE have also been simulated using the GCE model (Li et ah 1999; Tao 
2003). Sui and Yanai (1986) provided the GATE large-scale forcing for the GCE model. 
The environmental conditions for SCSMEX, TOGA COARE and GATE can be found in Tao 
et al. (2004). A TOGA-COARE surface flux algorithm (Wang et al. 1996) is used to 
calculate sea surface fluxes for these oceanic cases. 


(b) Continental Cases (ARM 1997 and 2002) 
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The Atmospheric Radiation Measurement (ARM) program set up the Southern Great Plains 
(SGP) site to observe clouds and precipitation for climate research. The site was centered at 
36.6°N and 96.5°W. Two field campaigns at the site, conducted in summers of 1997 and 
2002, are referred to here as ARM-SGP-97 and 02. The ARM forcing data were also 
obtained from the variational analysis approach of Zhang and Lin (1997) and Zhang et al. 
(2001). The surface fluxes are obtained from site-wide averages of observed fluxes from the 
ARM Energy Balance Bowen Ratio (EBBR) stations and are used to drive the model (Zeng 
et al. 2007). The ARM-SGP-97 numerical simulation starts at 2330 UTC 18 June 1997 and 
lasts for 29 days. The ARM-SGP-02 one starts at 2030 UTC 25 May 2002 and lasts for 20 
days. For the ARM cases, the surface fluxes are prescribed, and the surface wind does not 
interact with the boundary layer. 

All SCSMEX, TOGA COARE, GATE, ARM cases were simulated by GCE model 
previously. Please refer to Das et al. (1999). Li et al. (1999) and Zeng et al. (2009), for 
GATE, Johnson et al. (2002) and Zeng et al. (2009) for the TOGA COARE cases, Tao et al. 
(2003b) and Zeng et al. (2008) for the SCSMEX cases and Xu et al. (2002), Zeng et al. 
(2007, 2009) for the ARM cases with regard to the temporal variation of the wind shear. 
Table 2 su mm arizes all GCE model integrations for these cases. 

2.3 TRMMPR 


The TRMM Precipitation Radar (PR) operates at a frequency 13.8 GHz. This is a 
moderately attenuating frequency so that one of the major objectives of the PR algorithms is 
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to estimate and correct for the path attenuation. Because the PR is a single-wavelength, 
single-polarization, non-Doppler radar, there are only a few methods available to estimate 
rainfall. These include the Hitschfeld-Bordan (HB) method, the surface reference technique 
or SRT, a "hybrid" of these, and the mirror-image technique. The rain estimation algorithm 
uses a hybrid of the Hitschfeld-Bordan method and the SRT to correct for attenuation and 
derive an estimate of the range profile of the radar reflectivity factor, dBZ. The rain rate 
profile is then calculated from the dBZ profile and an appropriate Z-R relationship. The 
algorithm also includes surface clutter rejection and an attempt to correct for effects of non- 
uniform beam filling. 

Convective-stratiform classification of the rain is based on a combination of the 
vertical and horizontal structure of the radar reflectivity field. The vertical profile is checked 
for the presence of a bright-band (melting layer) by considering the behavior of the second 
derivative of the radar range profile. When the second derivative exceeds a threshold, a 
bright band is taken to be present and the rain is classified as stratiform. In cases where a 
clearly defined melting layer is absent, the horizontal rain structure is examined by means of 
a modified version of an algorithm designed for the analysis of ground-based radar data. 

Statistics of the instantaneous, high-resolution rain rates are compiled on a monthly 

basis over 0.5° x 0.5° latitude-longitude grids. Near-surface rain rates and the rain rates at 3 
altitude levels are stored according to rain type. The statistics include means, standard 
deviations and histograms of the rain rate, radar reflectivity factors, bright-band and storm 
heights. 
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3. Results 

3.1 PR estimated relationship between rainfall intensity and stratiform percentage 

Figure 1 shows the relationship between rainfall intensity and stratiform precipitation 
percentage for ocean and land region from PR. The difference between two regions is not 
significant. For example, there are more stratiform rain amount (percentage) for light rain 
(intensity) from 10 to 120 mm/day (or less than 1 mm/hour to 5 imn/hour) for both ocean and 
land region, but more convective rain amount for heavy rainfall (over 200 mm/day). This 
result might imply that there is general relationship between rainfall intensity and startiform 
amount. It also suggests that convective - stratiform heating algorithm could use this 
relationship to build its look-up heating table (more discussion in section 4). However, there 
are differences between oceanic and land region that showed in the low stratiform (or liigh 
convective) regime. For example, more convective rain for precipitation rate over 100 
mm/day for the land region than its oceanic counterpart. 

Figure 2 shows TRMM PR observed relationship between rainfall intensity and 
stratiform precipitation percentage over SCSMEX / NESA for May - June of 1998. It also 
shows the same SCSMEX region and the same period but for 10 years of observation. 
Again, the light (heavy) rainfall regime shows more (less) stratiform rainfall amount. Since 
SCSMEX / NESA region is where monsoon occurs and that may explain the more 
convective rain amount for heavy rainfall (similar to land as shown in Fig. 1(a)). Table 3 
shows PR observed rainfall amount and stratiform percentage for global tropics (52.8%) and 
its oceanic (55.0%) and continental (47,4%) components. More stratiform amount (about 7% 
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difference) was observed for oceanic than continental convection. Similar result also found 
in regional (SCSMEX-monsoon and ARM region). More stratiform rain amount occurs for 
SCSMEX than ARM cases. 

Schumacher and Houze (2003) examined the stratiform rain in the tropics using 
TRMM PR. Their results showed that stratiform accounts for 40% of total rain amount over 
20° N to 20° S for 3-year period (1998-2000). Their results also showed that convective rain 
rates attain greater value over land (from 26 to 38%) than ocean (from 41 to 53%) regions. 
In addition, their results found that 53% (36%) of stratiform rain amount occurs over 
Southeast Pacific (South and North America). Generally, these conclusions are good 
agreement with our results that are based on 10-year period of TRMM PR and over larger 
geographic area, except that more global tropic stratiform amount was shown in our results. 

3.2 GCE model simulated relationship between rainfall intensity and stratiform 
percentage 

In order to compare the TRMM 0.5 x 0.5 degree rainfall products, the GCE model simulated 
results were divided the domain (512 km) equally into 8 sub-domains so that each has a 
horizontal scale of 64 km. Then, average rainfall rates/amounts in convective and stratiform 
regions (including its percentage to total rainfall) in each sub-domain were calculated. Using 
the stratiform percentage and the average rate, the sub-domain is classified to be an element 
of a 36x20-size array, where the array has a precipitation rate interval of 20 mm/day (about 1 
mm/hour) and a stratiform percentage interval of 5%. Finally, calculate the vertical 
distributions of heating in each sub-domain and over the entire modeling period. 
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Figure 3 shows the model simulated PDF of the relationship between rainfall 
intensity and startifonn percentage over land and ocean regions. For oceanic region, a total 
of 123 day model integration combining SCSMEX, TOGA COARE and GATE cases. 
Shorter model integration for ARM cases (a total of 49 days by combining 1997 and 2002) 
was for land region. That is why a smother PDF was shown in the oceanic case. The model 
results also show similar feature as in PR observation (Fig. 1). For example, there is more 
stratiform rain amount (percentage) for light rain (intensity) from 10 to 120 mm/day (or less 
than 1 mm/hour to 5 mm/hour) for both ocean and land region, but more convective rain 
amount for heavy rainfall (over 200 mm/day) for land case. However, there is one major 
difference between PR observation and model simulation that is model simulates high 
stratiform rain amount for very small rainfall case (10 mm/day or < 1 mm/hour). The 
sensitivity of PR for not detecting very light rain (17 dBz) is the main reason for this 
difference. 

Figure 4 shows the relationship between rainfall intensity and startifonn percentage 
for each individual land and ocean case. Clearly, the two summer land cases (or ARM 1997 
and ARM 2002) show similar behavior in tenns of more convective rain for high rainfall 
intensity just as the oceanic cases. In addition, the land cases have higher frequency in heavy 
rainfall rate and small stratifonn percentage than the oceanic ones, which is consistent with 
the observations shown in Fig. 1 . 

Table 4 shows grid-averaged total rainfall and stratifonn percentage for GCE 
simulated cases. More rainfall is occuned in the oceanic cases than those of continental 
cases. This is due to the fact that more precipitable water (or moist) in oceanic environment 
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than its continental counterpart (see Table 1 in Tao et al 2004). The vertically integrated 
water vapor contents are very moist for the SCSMEX cases (over 62 g cm' - ) compared to the 
TOGA COARE and GATE cases. That is the reason for largest rainfall for SCSMEX 
simulation. For TOGA COARE and GATE cases, the environmental condition is more 
moisten in TOGA COARE but with less rainfall amount than GATE. It is because that the 
model simulation starts in November that did not have many active convective events. 
Generally speaking, more stratiform percentage (40-50%) should be simulated in the tropical 
oceanic cases than mid latitude continental cases. However, the ARM cases also have large 
stratiform rain amount/percentage (from 36 to 41%). Houze (1977), Zipser et al. (1981) and 
Gamache and Houze (1983) estimated that the widespread stratiform rain accounted for 
about 32%-49% of the total rainfall from the GATE cases. The fraction of stratiform rainfall 
from midlatitude squall lines has been estimated at 29%-43% (Rutledge and Houze 1987; 
Johnson and Hamilton 1988). The GCE model simulated results are in good agreement with 
these observations. 

By comparing GCE model simulation and PR observations (Table 2), the simulated 
rainfall amount for ARM (2002) case is much larger than the PR estimated. The sampling 
may be the reason for this difference. On the other hand, the GCE model simulated 
stratiform percentage for SCSMEX and ARM (2002) cases are in excellent agreement with 
TRMM PR estimated. In addition, the GCE model simulated rainfall amounts for SCSMEX, 
ARM (1997 and 2002) and GATE cases are in good agreement with those estimated by 
sounding networks. The TOGA COARE case is 17% less than sounding estimated. The 
GCE model has simulated two TOGA COARE cases (10-17 December 1992 and 19-27 
December 1992) and the simulated rainfall is an excellent agreement with sounding 



estimated (see Table 5 in Tao et al. 2004). We believe that the difference mainly occurs 
during the November 1992 model integration. 
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3.3 GCE model simulated each component of heating budget 

Qj can directly relate to the contributions of cloud effects which can be explicitly estimated 

by GCE model (Soong and Tao 1980; Tao and Soong 1986; Tao et al. 1993a and many 
others): 

& = 7r[- ( p " ° - V • VO'] + — [Lv(c -e) + L f (f - m)+L*{d - j)] + Qr (1) 

pdz c P 

where the primes indicate deviations from the large-scale enviromnent mainly due to small 
scale cloud processes. The variable 6 is potential temperature; p is air density; 
R/ c 

n — (p / P OQ ) p is non-dimensional pressure (where p and p OQ are dimensional and 
reference pressures with p QO taken as 1 000 hPa), and c p and R represent the specific heat at 
constant pressure and the gas constant of dry air, respectively. The variables L v , Lf and L s 
are the latent heats of condensation, freezing, and sublimation, respectively, while the 
variables (c, e , f m, d, s) denote rates for condensation of cloud droplets, evaporation of 
cloud droplets and rain drops, freezing of water droplets and rain drops, melting of ice 
crystals, snow flakes, graupel and hail, deposition of ice crystals, and sublimation of all ice 
hydrometeors, respectively. 

The term (1 / c p ) [^(c-e) +Lf(f -m) +L s (d -s)] is the LH due to microphysical 
phase changes. As defined in Yanai et al. (1973), Qj is the apparent heat source, while Q R is 



the radiative heatiiig/cooling rate associated with radiative transfer processes. The first two 
terms oil the right-hand side of Eq. (1) are the vertical and horizontal eddy heat flux 
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convergence ( 71 [ dw'6' / dz ] and 7r[V • V' 6' ] ), where the horizontal eddy term is neglected 
when Eq. (1) is spatially averaged over a large area suitable for large-scale diagnostic 
analysis or over the GCE model if a cyclic lateral boundary condition is applied. The 
horizontal eddy term (entrainment and detrainment) could not be negligible over a small area 
or over the sub-domain averaged in this study. 

Figures 5 and 6 illustrate each component of simulated heating structures associated 
with land and oceanic cases, respectively. These figures also show the characteristics of each 
heating component at specific stratifiorm amount (less than 5%) with different rainfall 
intensity (less than 20, between 80-100, 180-200, 380-400 and 480-500 mm/day) and at 
specific rainfall intensity (140-160 mm/day) with different stratiform amounts (less than 5, 
between 20-25, 45-50, 70-75 and 95-100%). These heating profiles are calculated over sub- 
domain (64 grids with 1 km spacing). Evident in the figures are: (1) latent heating is larger 
with large surface rainfall (Figs. 5a and 6a), (2) the level of heating maximum is higher with 
larger rainfall intensity (Figs. 5a and 6a), (3) a low level heating is evident for light rainfall 
(associated with shallow convective convection, Figs. 5a and 6a), (4) both eddy and radiative 
heating are quite smaller compared to heating but they cannot be negligible with light surface 
rainfall intensity (i.e., less than 100 mm/day, Figs. 5b, 5c, 5e, 5f, 6b, 6c, 6e and 6f), (5) the 
level of heating maximum is higher with less stratiform amount for oceanic case (Fig. 5d), 
(6) the cooling and heating aloft occurs at low troposphere with large stratiform amount with 
moderate rainfall (Figs, 5c and 6c). There are a few differences between land and ocean 
cases. For example, there are two heating maxima at low rainfall and small stratiform 
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amount for land case (Fig. 6a). The level of maximum heating is also higher for land case 
than oceanic case (Figs. 5a and 6a). In addition, the heating and eddy term is stronger for 
land case than oceanic case. This is because the convective vertical velocity is stronger for 
land case. 

These GCE model simulated heating profiles are in many ways in good agreement 
with observed convective and stratiform heating structures (Johnson 1984; Houze 1989, 
1997). These modeled convective and stratiform heating profiles will be in the new CSH 
algorithm's look-up table. 

4. Performance of consistency check of new convective and stratiform heating 

algorithm 

Current CSH algorithm uses surface precipitation rates and amount of stratiform rain. The 
CSH algorithm also utilizes a lookup table that consists of convective and stratiform diabatic 
heating profiles for various types of cloud systems in different geographic locations. These 
profiles stored in the lookup table are obtained from GCE model simulations or from budget 
calculations by temporally and spatially averaging the heating distributions in the convective 
and stratiform regions of the systems, which are then normalized by their total surface 
rainfall. Only about 15 heating profiles are in the current CSH heating algorithm. The 
improved CSH algorithm will also use the surface precipitation rates and amount of 
stratiform rain. The only difference is that much more heating profiles will be stored in the 
look-up table for improved CSH algorithm. These profiles will be selected based on the 
rainfall intensity and its stratiform amount. Some of heating profiles are shown in Figs 5 and 
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6 as discussed in Section 3.3 (only partially shown). In this section, the performance 
(including consistent check and comparison with previous CSH algorithm) of the improved 
CSH algorithm will be discussed. 


4.1 Consistency Check 

Consistency checking involving CRM-generated heating profiles and both algorithm- 
reconstructed and diagnostically-estimated heating profiles is a useful step in evaluating the 
performance of a given LH profile algorithm. In this process, the CRM simulation of a time- 
dependent precipitation process (multiple-day time series) is used to obtain the required input 
parameters for a given LH algorithm. The algorithm can then reconstruct the heating profiles 
that the CRM simulation originally produced, and finally both sets of conformal estimates 
(model and algorithm) can be compared to coincident estimates of diagnostically-based 
heating derived from radiosonde observations. Such observations from various field 
experiments, as well as simulations of individual precipitation systems, have been used for 
such consistency checks (Tao et al. 1990, 1993, 2000; Olson et al. 1999, 2006; Sliige et al. 
2004, 2006). 

Figure 7 shows the time series of the GCE model simulated, new CSH reconstructed 
heating and sounding estimated heating for the SCSMEX case. For the CSH reconstructed 
heating is using the GCE model simulated surface rainfall and its stratiform percentage. 
Also only the heating associated with raining regions (or boxes) are retrieved. The results 
show that the pattern in temporal variations of the SLH algorithm-reconstructed heating 
profile agrees well with that of the GCE-simulated heating profile. The GCE model 



simulated heating is stronger than sounding estimated. This is due to the fact that GCE 
simulated heating is based on five-minute cloud statistics (latent heating, eddy transport and 
radiative heating/cooling) compared to 3-hourly sounding estimates. 
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There are several differences between model simulated, CSH retrieved and sounding 
estimates. The CSH retrieved heating is generally slightly weaker than both simulated and 
sounding estimated. In addition, the retrieved heating structure does not show cooling due to 
the fact that the retrieved heating is only for raining region. There is no retrieved heating for 
those periods without surface rain (the radiative heating is much weaker than latent heating, 
see Fig. 5). Another major difference is the retrieved heating structure shows light heating 
above 10 km level (Steve and Xiping: Can we have explanation). 

There are strong heating and cooling at upper troposphere simulated by the GCE 
model. It is mainly caused by the radiative processes. For example, a cloud emits infrared 
radiation at top upwards and absorbs the upward infrared radiation at bottom, which results 
in radiative cooling and heating at top and bottom, respectively. Since the heating at bottom 
increases instability and therefore large eddies there, the heating region is extended upwards. 
In contrast, the cooling at top increases instability and large eddies right below the cooling 
region, and thus the heating from eddies offsets the radiative cooling partly. 

Forty four-day-averaged profile of the GCE model simulated, the CSH reconstructed 
and sounding estimated Q1 are shown in Fig. 8. These profiles are simply the time averages 
of the profiles shown in Fig. 7. The 44-day-averaged profiles of the reconstructed heating 
agree very well with the GCE simulated although the reconstructed heating is a little bit 
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weaker than the simulated one. The level of maximum heating simulated and reconstructed 
is about 7 km altitudes and is in good agreement with sounding estimated. Both simulated 
and reconstructed heating profiles have a distinct cooling near 4. km, which is due to the 
melting processes. This is because the eddy heat flux convergence compensates for the 
distinct cooling due to the melting. This feature is also shown in SLH reconstructed (see 
Sliige et al. 2004). Both GCE simulated and CSH reconstructed is stronger than sounding 
estimate. Note that the model did not simulate strong cooling in middle and low troposphere 
between active convective events (Fig. 7). 

4.2 Performance 

The upper two panels of the Fig. 9 illustrate 10-year mean apparent heating ( Q ] ) at two 
different altitudes (0.5, 4 and 7 km) over the global Tropics from the new CSH algorithm 
based upon the PR-monthly rainfall product. The new CSH algorithm has new heating look- 
up table that is binned based on rainfall intensity (mm/day) and stratiform amount 
(percentage). Then, the Q ] profiles are calculated based on PR-retrieved surface rain rates 
and stratiform fractions at a grid scale of 0.5 degrees for either oceanic or continental 
locations. The concept between the new and old CSH algorithm is same. The only major 
difference between the new and old CSH algorithm is heating profile stored in the look-up 
table. The old CSH algorithm only consists of normalized convective/stratifonn kernel 
profiles are created by averaging all 16 (4) oceanic (continental) base profiles in the CSH 
algorithm’s complete library of 20 heating profile pairs distributed regionally and according 
to storm type. The new CSH algorithm consists more than 200 heating profiles in the 
oceanic and continental region, respectively. 
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As expected from the design of the CSH algorithm, the horizontal distribution of the 
estimated Qj structure is similar to the pattern of surface rainfall (lower panel of the Fig. 8). 
For example, well defined ITCZs in the east and central Pacific and Atlantic Oceans, a well- 
defined South Pacific Convergence Zone (SPCZ) in the central-southern Pacific Ocean, and 
broad areas of precipitation events spread over the continental regions are all evident. Also, 
strong LH release in the middle troposphere (5 °C day" 1 and greater) is associated with 
heavier surface precipitation. Heating in the upper troposphere is weaker than those in the 
middle troposphere. 

Figures 10(a) and 10(b) are the 10- year averaged Q1 derived from old and new CSH 
algorithm. Clearly, there are two heating maxima derived from both old and new CSH 
algorithm. Both occur about the same locations, and respond to the surface rainfall 
associated with the ITCZ (at 5° North) in the North Hemisphere, and one associated with the 
SPCZ (2.5° Soudi). The heating is stronger in the North Hemisphere because the surface 
rainfall is larger in the North Hemisphere. This difference could have impact for large-scale 
circulation respond to the heating gradient between North and South Hemisphere. One major 
difference between old and new CSH algorithm is the new CSH retrieved stronger heating 
than that old CSH. The level of maximum of the heating is lower in the new CSH algorithm 
derived. Another major difference between old and new CSH algorithm is at the low 
troposphere. The heating at low troposphere is stronger in the new than old CSH algorithm 
near the equator. 
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Figure 10(c) shows the mean 10-year old and new CSH-retrieved average heating 
profile over the Tropics. The land and ocean components are also shown from both new and 
old CSH algorithm for comparison. There are very little variation of heating between land 
and ocean for old CSH CSH algorithm. For new CSH algorithm, the heating is much 
stronger between 2 and 8 km altitude over the land than over the ocean. The stronger 
convective heating in the land could be the consequences of stronger vertical velocity and / 
or large convective rain amount (Table 3). In addition, there is cooling at low troposphere in 
the new CSH over land. The drier environment in the land than the ocean could allow more 
evaporation cooling over the land region. Another major difference between new and old 
CSH algorithm is that level of maximum heating is lower (about 6 km altitude) for the new 
CSH algorithm. 

Lin and Johnson (1996) compared the mean latent heating from the TOGA COARE, 
Marshall Islands region (Yanai et al. 1973) and the GATE region (Thompson et al. 1979). 
Their results show that the level of maximum heating is about 6-6.5 km altitude for the 
Western Pacific (TOGA COARE and Marshall Island region), but is only about 4 km altitude 
for the GATE region. The lower level of maximum heating the GATE region may be due to 
lower SST in eastern Atlantic (Thompson et al. 1979). Greco et al. (1994) calculated latent 
heating profiles over S. America from the sounding network network. Their results indicated 
that the distribution of heating is quite similar to the studies of those of West African squall 
lines (Chong and Hauser 1990). Peak heating also occurs between 500 and 550 hPa (about 
5-6 km). The new CSH algorithm retrieved level of the maximum heating is also about 6 km 
altitude that is in good agreement with these diagnostic studies. However, the level of the 
maximum heating is higher (about 7.5 km) for monsoon over S. China Sea (SCSMEX 
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NESA, Johnson and Ciesielski, 2002) and for convective system over S. America (LB A, 
Halverson et al. 2002). The regional analyses of the new CSH products including 
comparison with diagnostic studies and with other heating products (i.e., products derived 
from other heating algorithms, SLH and TRAIN) will be conducted in next paper. 

5. Summary 

The relationship between rainfall intensity and stratiform cloud was examined by analyzing 
TRMM PR and GCE model simulated results. The relationship between surface rainfall 
intensity, stratiform percentage and latent heating was also examined by using GCE model 
simulated precipitation processes. These results were used for improving the CSH 
algorithm. The performance of the new CSH algorithm was presented by comparing the 
results with previous CSH algorithm. The major highlights are as follows: 

. Both PR estimated and GCE model simulated results showed that are more stratiform 

rain amount for light rain for both ocean and land region, but more convective rain amount 
for heavy rainfall. The differences between oceanic and land region is that more convective 
rain for precipitation rate over 100 mm/day for the land region than its oceanic counterpart. 

. One major difference between PR and GCE simulated rain intensity and stratiform 

relationship is in very light (less than 10 mm/day) rainfall region. This difference is caused 
by the PR sensitivity (17 dBz) on light rainfall. 

. The GCE model simulated the stratiform percentage is in excellent agreement with PR 


estimated for SCSMEX (1998) and ARM (2002) cases. But the GCE simulated more rainfall 



amount for ARM case compared to PR estimated. On the other hand, the GCE model 
simulated rainfall amount is in good agreement with sounding estimated. Sampling is the 
main cause for the difference between PR estimated and GCE model simulated rainfall 


24 


amount. 

. Each component and profile (latent heat, eddy term and radiative/cooling/heating) of Q1 

associated different stratiform percentage and rainfall intensity for land and ocean cases, was 
simulated by GCE model. These profiles are used for new CSH algorithm’s look-up tables. 
There are many similarities in these profiles between land ocean cases. For example: (1) 
latent heating is larger with large surface rainfall, (2) the level of heating maximum is higher 
with larger rainfall intensity, (3) a low level heating is evident for light rainfall, (4) both eddy 
and radiative heating/cooling are quite smaller compared to heating but they cannot be 
negligible with light surface rainfall intensity, and (5) the cooling and heating aloft occurs at 
low troposphere with large stratiform amount with moderate rainfall. The major differences 
are that there are two heating maxima at low rainfall and small stratiform amount for land 
case and the level of maximum heating is also higher for land case than oceanic case. In 
addition, the heating and eddy term is stronger for land case than oceanic case. 

. One major difference between new and old CSH algorithm is that the level of 

maximum heating is 1 to 1.5 km lower for the new than old CSH algorithm. Another major 
difference is that there is more variation in heating profiles in the ocean and land regions. For 
old CSH algorithm retrieved heating: there is almost no variation of heating between land 
and ocean for old CSH algorithm. 
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The results from this paper are two-dimensional version of GCE model was used to 
provide the heating profiles for CSH algorithm. It is known that the results from 2D and 3D 
CRM simulations can differ not only in cloud dynamics but also microphysics (i.e., Tao et al. 
1986; Johnson et al. 2002; and Zeng et al. 2008). More convective cores with stronger 
vertical velocity exist in a 3D CRM than in a 2D one. As a result, the cloud ice content in 
the upper troposphere in a 3D CRM is higher than that in its 2D counterpart (Zeng et al. 
2008). The sensitivity of the relationship between rainfall intensity and its stratiform amount 
using 3D needs to be examined. 

Only limited CRM simulated cases were used in this study. The Observations from 
additional field experiments (TWP-ICE, GAME and future GPM Ground Validation (GV) 
site) will be needed to provide new types of initial conditions for GCE model. In addition, 
the large scale analyses (i.e. NASA Goddard Modern era retrospective-analysis for research 
and applications, MERRA) could also provide initial and condition and forcing for CRM. 
Figure 10 shows the preliminary results using MERRA to provide the forcing for GCD 
model. The results indicate that the GCE+MERRA approach can provide reasonably good 
simulations of reality. The two simulations generally agree with each other to a large extent. 
Thus, the GCE+MERRA, are compatible and could be used for additional cases (especially 
for those regions with large surface rainfall, such as Indian Ocean, SPCZ, S. America, and 
Africa) for CSH algorithm. 
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TRMM Data 
Needed 

Heating 

Products 

Key References in Algorithm 
description 

CSE[ (Convective- 
Stratiform Heating) 

PR, TMI, 
PR-TMI 

Qi, LH 

Tao et al. (1993; 2000; 2001) 

SLH (Spectral Latent 
Heating) 

PR 

LH, Qj-Qr 

Shige et al. (2004; 2007, 2009) 

TRAIN (Trained 
radiometer algorithm) 

TMI (PR 
Training) 

Qi-Qr, LH 

Grecu and Olson (2006), Olson et 
al. (2006), Grecu et al. (2009) 

HH 

(Hydrometeor Heating) 

PR-TMI 

LH 

Yang et al. (1999; 2006) 

PRH (Precipitation 
Radar Heating) 

PR 

LH 

S atoll and Noda (2001), 
Katsumata et al. (2009) 


Table 1 Summary of the five LH algorithms participating. Data inputs, retrieved 
products, and salient references included. Note conventional relationship 
between Qi (apparent heat source), LH, and Q R ( radiative heating) is expressed 
by Qi-Qr = LH + EHT, where the final term represents eddy heat transport by 
clouds (noting that vertically integrated EHT is zero, i.e., it provides no explicit 
influence on suiface rainfall). Note that CSH, SLH and TRAIN explicitly used the 
CRM simulated latent heating profiles as their heating algorithm - look up tables. 
Both HH and PRH also applied implicitly CRM simulated results (i.e., cloud 
vertical velocity). 


Field Campaign 

Geographic Location 

Starting Date 

Modeling Days 


ARM-SGP-97 

(37°N, 97° W) 

18 Jun 1997 

29 

Tao et al. (2004); 
Zeng et al. (2009) 

ARM-SGP-02 

25 May 2002 

20 

Zeng et al. (2007, 
2009) 

SCSMEX/NESA 

(21°N, 1 1 7°E) 

6 May 1998 

44 

Tao etal. (2003), 
Zeng et al. (2008) 

TOGA-COARE 

(2°S, 1 54°E) 

1 Nov 1992 

61 

Das et al. (1999); 
Johnson et al. (2002); 
Zeng et al. (2009) 

GATE 

(9°N, 24°W) 

1 Sep 1974 

18 

Li et al. (200x); 
Zeng et al. (2009) 


Table 2 Field campaigns (ARM, SCSMEX, TOGA CO ARE and GATE) including the 
geographic location, starting time and integration length of GCE model simulations. 
Also include the previous GCE modeling papers that have simulated the case. 
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Rainfall amount (mm/day) 

Stratiform Percentage (%) 

Global Tropics 

2.26 

52.8 

Ocean 

2.35 

55.0 

Land 

2.06 

47.4 

SCSMEX (1998) 

11.24 

42.3 

SCSMEX (10 years) 

3.64 

49.8 

ARM (2002) 

2.59 

35.7 

ARM (10 years) 

2.36 

41.1 


Table 3 PR estimated rainfall amount and stratiform % (Global, Land, Oceans, SCSMEX 
(1998), SCSMEX (10 year), ARM (1997, 2000) and ARM (10 year)) 



Simulated Rainfall 
amount (mm/day) 

Stratiform 
Percentage (%) 

Estimated Rainfall 
amount (mm/day) 

SCSMEX 

12.31 

42.6 

11.35 

ARM (1997) 

4.31 

41.3 

4.32 

ARM (2002) 

4.85 

36.0 

4.77 

TOGA COARE (1992-1993) 

7.72 

47.6 

9.32 

GATE (1974) 

10.56 

41.4 

11.38 


Table 4 CRM simulated rainfall amount and stratiform % (SCSMEX (1998), ARM (1997, 
2002), TOGA COARE (1992) and GATE (1974)) 





Stratilorm Precipitation Percentage (%) 
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Fig. 1 Hit frequency of the TRMM observations over (a) oceanic and (b) continental 
regions. 



Fig. 2 


(a) The same as Fig. 1 except for TRMM observations over SCSMEX/NESA for 
May- June of 1998. (b) as (a) except for 10 years of observations. 


Stratiform Precipitation Percentage [%) 


35 


Model Ocean Model Land 
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Fig. 3 Same as Fig. 2 except for GCE model simulation, (a) is for oceanic (TOGA 
COARE, GATE and SCSMEX) case and (b) is for continental (ARM 1997 and 
2002) case. 
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Model TOGA Model ARM-SGP-97 



Model GATE 



Precipitation Ra le (mm/day) 
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ARM-SGP-02 


Model SCSMEX/NESA 



Fig. 4 Same as Fig. 2 except for GCE model simulation. Left three panels are oceanic 
cases, TOGA COARE, GATE and SCSMEX. Right two panels are continental 
cases, ARM 1997 and 2002. 
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CSH LUT LH - OCEAN 



CSH LUT LH - OCEAN 



CSH LUT Qla - OCEAN 


CSH LUT Qla - OCEAN 




CSH LUT Qlr - OCEAN CSH LUT Qlr - OCEAN 




5 Vertical profiles of the latent heating rate (top), eddy rate (middle) and radiative 
rate (bottom). Line labels indicate the range of the convective precipitation 
percentage. 
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Q1 (K/day) Model 
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Fig. 7 Time-height cross sections of the modeled (top), retrieved (middle) and observed 
(bottom) Q1 in SCSMEX/NESA case. 
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Fig. 8 


Vertical profiles of mean Q1 over SCSMEX/NESA that is observed (thick 
solid), modeled (thin solid), retrieved (thick dashed line). 
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Fig. 9 Ten-year mean Qj heating rates at 7, 3, and 1 km AGL (upper 3 panels) along with 
surface rain rates (lower panel) over the global Tropics determined by the Goddard 
Space Flight Center (GSFC) Convective-Stratiform Heating (CSH) algorithm 
applied to 1998-2002 Precipitation Radar (PR) measurements acquired from the 
Tropical Rainfall Measuring Mission (TRMM) satellite 
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Fig. 10 Zonal distribution of heating structure from new (top) and old (middle) CSH 
algorithm. The bottom shows the global averaged heating profiles (including land 
and ocean component using new and old CSH algorithm). 
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Fig. 11 The bold black line is the 3-hourly time series of domain-averaged precipitation 
observed during SCSMEX period (from 0600 UTC 18 May to 0600 UTC 26 May 
1998). The red line is the 3-hourly precipitation time series from GCE simulation 
driven by SCSMEX intensively observed large-scale tendencies. The green line is 
same to the red line except the GCE was driven by tendencies from MERRA 
output. The two simulations generally agree with each other to a large extent and it 
indicates that the GCE+MERRA approach can provide reasonably good simulations 
of reality. 


